The issue of aliasing in PSR B0943+10 
II Signal processing arguments 

M. Vivekanand* 
National Center for Radio Astrophysics, TIFR, 
Pune University Campus, P. O. Box 3, 
Ganeshkhind, Pune 411007, India. 

February 5, 2008 



Abstract: DeshDande & Rankin (1999, 2001) claim that the frequency of the very 



narrow feature, in the spectrum of radio flux variations of PSR B0943+10, is an 
alias of its actual value. They also claim to have detected an amplitude modulation 
on the above phase modulation. This paper argues that both these claims are 
unjustified. 
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1 Introduction 

The rotation powered radio pulsar PSR B0943+10 exhibits sub pulses that are 
drifting systematically from period to period within the observable pulse window. 
This pulsar has a very narrow feature in the longitude resolved spectrum of in- 
tensity fluctuation s; its Q, defined as its cen t ral frequency divided by its width 



is relatively high (Taylor & Hueguenin 1971; Backer. Rankin & Campbell 1975; 



Sieber & Osteiill975h . Recently iDeshpande & Rankinl (1999. 2001; henceforth 
DR1999 and DR2001) put its Q at > 500. It has been speculated that this very 
narrow spectral feature, occurring at 0.465 cycles per pulsar period [cpp], could 
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be an alias of the actual value 0.535 [cpp] (iSieber & Ostenll975l) . Backed d 1973b 
noted such a general possibility in radio pulsars and claimed, based on the phase 
information in the fluctuation spec trum, t hat PSR B2303 +30 has an aliased spec- 
tral feature. This was contested bv lSieber & Osten (119751) . who s tate on their page 
326 that "even a phase analysis, contrary to Backer's statements dBackeitil973L is 
unable to decide between the two possibilities", viz., whether the spectral feature 
is aliased or not. Indeed, in a later paper Backer did not repeat such a claim for 
PSR B0943+10, whose fluctuation spectrum is similar to that of PSR B2303+30 
(Bac ker. Rankin & Campbellll 19751) : he quoted the the true value as 0.465 [cpp]. 
However, DR1999 and DR2001 recently claimed that the spectr al feature is in- 
deed a n alias; that the actual value is 0.535 [cpp]. In their view ISieber & O ster 
( 1975) "came to the wrong conclusion" (DR2001). They also claim that the weak, 
symmetrically spaced sidebands, at 0.027 [cpp] away from the above spectral fea- 
ture "strongly suggest ... a regular, highly periodic amplitude modulation of the ... 
drifting sub pulse sequences". 

This paper argues that DR1999 and DR2001 are unjustified in drawing these 
two conclusions. The question, whether PSR B0943+10 has an aliased spectral 
feature or not, is still unresolved; and the latter observation is as likely, if not more 
likely, due to an additional phase modulation of the drifting sub pulses. 



2 A review of the relevant signal processing 



The signal from an ideal pulsar with drifting sub pulses falls under the topic "pulse 
position modulation" (PPM); it consists of periodically occurring narrow pulses 
whose positions are modulated by another periodicity. Its general p rinciples ca n 
be found in books on electronic communication engineering (see lLathil 1998). 
The spectrum of a PPM signal when t he posit i on mo dulation is due to a pure 
"tone" (a single frequency) is g iven bv IPanten (119651) on his page 541 and by 
Schwartz . Benn ett & Steirl d 19661 : henceforth SBS1966) on their pages 252 - 253. 



2.1 Frequency domain discussion 

To begin with let us assume that the drifting sub pulse pattern is PPM due to a 
pure tone. 

The abscissa in fig. 6-2-3 on page 251 of SBS1966 represents the phase of the 
sampling signal at which a pulse is observed; in our case the sampling signal has 
the pulsar period P (sec), or frequency 1 /P Hz or 1 [cpp]. The ordinate represents 
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the corresponding phase of the modulating signal; in our case it has period P3 
pulsar periods (frequency 1/P3 [cpp]) which is the repetition time of the drifting 
sub pulse pattern. The motion in time of the drifting sub pulses in this figure is 
along a straight line of slope 1 /P3. In the last para on their page 254, SBS1966 
state "the average pulse repetition frequency should be at least twice the highest 
signal frequency in order to obtain the minimum number of samples necessary for 
satisfactory signal recovery"; i.e., 1 > 2/P3 (=> 1/P3 < 0.5) to avoid aliasing. 
This can be verified by considering two straight lines of slopes < 0.5 and > 0.5 in 
fig. 6-2-3 of SBS 1966. 

Thus, the Nyquist sampling criterion is the same for a PPM signal and a canon- 
ical pulse amplitude modulation (PAM) signal (amplitude modulation of periodic 
pulses). DR2001 are wrong when they claim in their conclusion that "A harmonic 
resolved fluctuation spectrum uses the information within the finite width of the 
pulse to achieve a Nyquist frequency of 1 [cpp], showing clearly that the primary 
feature is aliased". Their average pulse repetition frequency is obviously P (sec); 
so their Nyquist frequency is only 0.5 [cpp]. Consequently, they are also wrong 
in concluding that "the primary feature is aliased" based merely upon the Fourier 
technique; they require additional and independent information, as discussed 
ahead. 

Consider a PAM signal in which the amplitude modulation is due to a pure tone 
of frequency v [cpp]. The amplitudes occur at frequencies v [cpp], 1 + v [cpp], 
2 + v [cpp], etc; and there is no difference in the amplitude spectra of the original 
(v < 0.5) and the aliased (v > 0.5) signal. However the phase spectra differ. 
Therefore one can distinguish between the original and aliased PAM signals only 
by using additional information such as the phase of the modulating signal. In the 
PPM case, the amplitudes in the spectra occur at frequencies v [cpp], l+mxv 
[cpp], 2+mXv [cpp], etc., where m = 1, 2, 3, ... is the order of the harmonic (see eq. 
6-2-13 on page 252 of SBS 1966). Now, both amplitude and phase spectra differ 
for the original and aliased signals. Therefore once again, one can distinguish 
between the original and aliased PPM signals only by using additional information 
such as the exact shape of the two amplitude spectra. In practice it is impossible 
to predict this exactly in the current pulsar context. 

Fig. H] simulates the amplitude spectra of a pulsar signal showing the drifting 
sub pulse phenomenon, modeled as a PPM signal with a saw-tooth modulation in 
time, for both the original modulating frequency (top panel), and its alias but with 
the opposite drift direction (bottom panel). The pulsar parameters are taken from 
DR2001 for PSR B0943+10; however the figure is insensitive to small variations 
in these parameters. The algorithm was checked by reproducing the spectrum in 
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Figure 1: Simulated amplitude spectra of a pulsar signal showing the drifting sub 
pulse phenomenon. Only a small range of frequency has been shown for better 
visual comparison. The time series consists of 4x 1024 x 1024 samples, each of du- 
ration 0.25 milli seconds (ms). The pulsar parameters are those of PSR B0943 + 10, 
taken from DR2001 : period P = 1 .0977 (sec), Gaussian integrated profile of width 
= 0.031 x P (sec), maximum time departure of drifting sub pulses = 0.021 x P 
(sec). The sub pulses are assumed to be Gaussian in shape of very narrow intrin- 
sic width, of w one sampling interval; making the width a more realistic value 
merely suppresses the spectra at higher frequencies. The drifting sub pulse phe- 
nomenon is modeled as a PPM signal with a saw-tooth modulation in time. Top 
panel: Original modulating frequency = 0.465 [cpp], or P3 = 2.1505376 periods; 
Bottom panel: aliased modulating frequency = 0.535 [cpp], or P3 = 1.8691589 
periods, and with opposite drift direction. 
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eq. 6-2-13 of SBS1966. The difference between the two amplitude spectra in 
fig- CD is generally < 5% of the maximum amplitude of pulsar harmonics, upto a 
frequency of 100 [cpp], for the harmonics m = 1 to 3. For example, at 59.93 
[cpp], the second harmonics differ by « 12.5% in the two panels; but their average 
amplitude is < 33% of the maximum amplitude. The third harmonics, aliased to 
59.605 [cpp], differ by w 9%; but their average amplitude is only about 4.5% of 
the maximum amplitude. At 110.395 [cpp], the harmonics differ by « 20%; but 
their average amplitude is now about 37.5% of the peak amplitude of the pulsar 
harmonics; this is also insignificant. 

Therefore at lower frequencies (< 100 [cpp]; the range analyzed by DR2001) 
the spectra of the original and the aliased signals should be predicted to an accu- 
racy of < 5% to tell the difference. This is almost impossible, considering other 
factors ignored here: random variations in shape, size and intensity from pulse 
to pulse, lack of knowledge of the exact modulating function, the exact shape of 
the integrated profile, etc. Maybe there is some hope of making better predic- 
tion about the higher harmonics (m > 3) at higher frequencies (> 100 [cpp]), but 
DR2001 do not analyze these frequencies. The next section indicates why, if at 
all, the two spectra can probably be distinguished only at higher frequencies. 

2.2 Time domain discussion 

Let t n be the time of occurrence of the n th pulse in a PPM signal that is modu- 
lated by a pure tone, of phase <p m , and frequency 1/P3 [cpp], implying angular 
frequency co m = 2n/(P3P), 

t„ + t sin (co m t n + <p m ) = nP, ( 1 ) 

where r is the maximum position departure (see eq. 6-2-2 of SBS1966). Let At n 
be the difference between t„ and nP, its time of occurrence in the un modulated 

case, 

At n = -t sin (O n + (O m At n ) , (2) 

where O n = 2nn/P3 + <p m . In PSR B0943+10, r « 0.021P, and the original 
P3 « 2.15, so the phase oj m T = 2ml(P3P) « 0.06, which is a small quantity. 
Furthermore, At n < r; therefore, 

At n « -t sin (O n ) [1 - r {0.5u) c - Aw,,,} cos (O,,)] (3) 

where a> c = 2n/P is the sampling frequency (angular) and Aco m = 0.5co c - a> m . In 
PSR B0943+10, Ato m (= 2.862 - 2.662 = 0.2) is much smaller than either 0.5o; c 
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Now consider a PPM signal that is modulated by a pure tone, of phase (p' m = 
-(p m , and a frequency that is the alias of the earlier one; i.e., of angular frequency 
(o' m = co c - co m . Let it have the same magnitude of drift but opposite in direction, 
i.e., t' = -r; Then it can be shown that 

At' n « -r sin (O n ) [1 + r {0.5 to c + AcoJ cos (O n )] , (4) 

for integer n. For non-integer n, which implies that the pulses start at some arbi- 
trary phase within the period, a constant phase adds to <D„ in eq. (4). 

Now, the observed time series corresponds to either eq. (3) or eq. (4). If one 
knows its absolute phase, then one directly compares the observation to those two 
equations, and thus determines if the observed frequency is original or aliased. 
However, this is rarely the case. So, one has to distinguish between the the 
time series +0.25u> c t 2 sin(20„) of eq. (3) and eq. (4). This is possible only if 
the sampling interval used for the observations is much smaller than the quantity 
0.5u) c x r 2 , which is w 1.52 ms in PSR B0943 + 10. However the sampling interval 
of DR2001 (1.006 ms) is barely smaller than this. Therefore it is probably impos- 
sible for DR1999 and DR2001 to distinguish between the original drift frequency 
and its alias with opposite drift direction. This argument also holds for saw-tooth 
modulation drifting pattern, which can be Fourier decomposed into harmonics of 
the fundamental modulating frequency a> m . 

Now one can understand the result of the previous section - a time resolution 
<< 1.52 ms implies a minimum Nyquist frequency of >> 0.5/0.00152 « 329 
Hz or 361 [cpp]. Indeed, the spectra in fig.[T]differ most in the range 500 to 900 
[cpp], mainly in the higher harmonics m » 3. For example, the 15 th harmonics 
(aliased) at frequency 679.975 [cpp] differ in amplitude by « 107%; and their 
average amplitude is about 24% of the maximum amplitude of pulsar harmonics. 

In summary, DR2001 can hope to distinguish between the original and aliased 
spectral feature in PSR B0943+10 by predicting the two amplitude spectra at fre- 
quencies much larger than 100 [cpp], particularly for the harmonics m > 3 to 15, 
and comparing those with observations. This they have not done. 

3 Aliasing claim by DR2001 

This section discusses the three main arguments that DR2001 appear to offer (their 
sections 3 and 5) regarding the spectral feature in PSR B0943+10. 
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3.1 Amplitude spectrum argument 

In the longitude resolved spectrum in fig. 1 of DR2001, the narrow spectral feature 
falls at frequency 0.465 [cpp]. The spectrum of such a signal will naturally have a 
Nyquist sampling limit of 0.5 [cpp], by definition. To test whether this feature is 
aliased or not, one needs data with much faster sampling. DR2001 achieve this by 
obtaining the spectrum of the original time series (including zero padding) which 
is sampled at 1.006 (ms). The corresponding spectrum certainly shows that the 
spectral feature is much stronger at 0.535 [cpp] than at 0.465 [cpp] (their fig. 4). 

Therefore they argue (on the right side of page 442) that "... the asymmetry 
between the 0.535 and 0.465 [cpp] features in fig. 4 is so great that we can regard 
them as a signature of an almost pure phase modulation". This is justified, as 
seen in fig. [T] above. Then, after a short paragraph offering no new argument or 
evidence, they state that "We can now be certain that the principal feature in fig. 1 
is the alias of a fluctuation the actual frequency which is greater than 0.5 [cpp] 
This is unjustified. 

It is purely fortuitous that the aliased feature is stronger than the original in 
their fig. 4. It depends upon the specific combination of true P3, true drift direc- 
tion and the exact shape of the IP of PSR B0943+10. Fig. El verifies the above 
argument. It shows the amplitude spectra of PSR B0031-07, another drifting pul- 
sar, observed using Ooty Radio Telescope (ORT). The top panel corresponds to 
fig. 3 of DR2001. In the bottom panel of fig. 0the data is altered: the period 
sequence is retained, but within each period the time sequence is reversed. This 
could have occurred either due to (1) PSR B0031-07 having the aliased drift fre- 
quency but the same drift direction, or (2) it having the original drift frequency 
but opposite drift direction. The asymmetry has also got reversed in the bottom 
panel of fig. El 

3.2 Phase spectrum argument 

DR2001 do not use the phase information corresponding to their fig. 4, but use that 
available in the longitude resolved spectrum (their fig. 1). They compare the rate 
of change of phase, with longitude, of the 0.465 [cpp] spectral feature, with that 
of two other much weaker spectral features occurring at 0.071 [cpp] and 0.607 
[cpp], which are supposed to be the aliased first and second harmonics of the 
fundamental. They claim on page 443 that "the harmonicity of these phase rates 
argues strongly that the 0.071 [cpp] and 0.607 [cpp] features are indeed the aliases 
of the second and third harmonics of the primary feature which is justified. 
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Figure 2: Top panel: Amplitude spectrum of 1024 periods of radio flux data 
of PSR B0031-07, obtained on 1996 March 31 using the ORT. The data was 
acquired in a gated fashion and later zero padded., as in DR2001. The sampling 
interval was 6.5 ms and the pulsar has period 0.9430 (sec). Only a small range of 
frequency has been shown for better visual comparison. Bottom panel: The same 
as in the top panel, except that the time sequence within each period is reversed, 
while the period sequence itself is retained. 
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However, this will be true irrespective of whether the original frequency is 
0.465 [cpp] or 0.535 [cpp]; this is obvious from Fourier theory. The only dif- 
ference is the inverted phase relation between the fundamental and its (aliased) 
harmonics, about which one has no independent information anyway. This author 
is confused about what DR2001 intended in invoking this discussion. This author 
is most confused about a sentence in the last para of their section 3, the section 
that is supposed to justify their claim of aliasing. It goes "The ambiguity, how- 
ever, between the remaining two possible combinations - namely, P3 < IP (thus 
implying negative drift) and P3 > IP (positive drift) ca n not be fully resolved 
through the above analysis, Presumably this is what ISieber & Oste r (1975) 
implied, with which DR2001 apparently disagreed earlier. It appears that either 
DR2001 are contradicting themselves, or this author has not understood what was 
it in PSR B0943+10 that was "otherwise multiply folded" which DR2001 have 
managed to "unfold" (last para of their section 3). 

3.3 "Modulation on modulation" argument 

DR1999 claim on their page 1010 that the symmetrical sidebands associated with 
the primary feature in the fluctuation spectrum represent "an amplitude modula- 
tion on the phase modulation". These sidebands fall 0.027 [cpp] higher and lower 
than the primary feature. Since 0.535/0.027 ~ 20, an integer value, DR2001 
claim in their conclusion that "it is at this point that we have conclusive evidence 
that the aliasing question is resolved 

Fig. |3] is similar to fig. [T] except that (1)4 times fewer samples and 4 times 
larger sampling interval, and (2) the drifting P3 = 3 periods, for better visibility. 
However, in the top panel the pulsar signal is additionally amplitude modulated 
by a sinusoid of frequency 0.05 [cpp]; of relative magnitude of modulation of 
50%, for the sidebands on the position modulating harmonics to show up signif- 
icantly. The pulsar harmonics also have the symmetric sidebands, but of much 
larger amplitude. In the bottom panel of fig.|3]the above signal is instead addition- 
ally position modulated, by a sinusoid of the the same frequency, and a maximum 
time departure of 0.010 [P]. Now also the pulsar harmonics have the symmetric 
sidebands, but much smaller in amplitude. 

From fig. |3] shows that both an additional amplitude modulation or an addi- 
tional position modulation of the pulsar signal in fig. [T] gives rise to sidebands 
around the drifting harmonics. However, in the former case the side bands are 
much stronger around the pulsar harmonics (at 1 [cpp], 2 [cpp], etc.), and in the 
latter case they are much stronger around the drifting harmonics, which is what 
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Figure 3: Same as in fig.[U except (1) 1024 x 1024 samples used, (2) sampling in- 
terval of 1.0 ms, (3) modulating P3 = 3 periods. Only a small range of frequency 
has been shown for better visual comparison. Top panel: The PPM signal is ad- 
ditionally amplitude modulated by a tone of period 20 [P], and relative amplitude 
0.5. Bottom panel: The PPM signal is additionally position modulated by a tone 
of period 20 [P], and maximum time departure of 0.010 [P]. 
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DR2001 observe. Therefore it is as likely, if not more likely, that DR2001 have 
noticed an additional position modulation over and above their drifting sub pulse 
pattern. DR2001 do not reason why they prefer the former over the latter. 

The foundation of DR2001 vanishes if they are unable to justify the amplitude 
modulation of the drifting sub pulses. 



4 Summary 

DR2001's claim of aliasing in PSR B0943+10 may rest purely upon an arbitrary 
assumption, that of a very steady "amplitude m odulation on the phase modula- 
tion". Recent claims ( Edwards & Stappersl 2002b to have independently verified 



the results of DR2001 are, in the opinion of this author, not relevant here because 
(1) they are model dependent, and (2) they work with the phase resolved spectrum, 
which is inadequate by definition to discuss the issue of aliasing. 

This research has made use of NASA's Astrophysics Data System (ADS) Bib- 
liographic Services. 
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